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We present results of quenched lattice calculations of the matrix elements relevant for Bd—Bd 
^ ' and Bs—Bs mixing in the standard model. Results for the corresponding S'C/(3)-breaking ratios, 

, which can be used to constrain or determine \Vtd\, are also given. The calculations are performed 

■ at two values of the lattice spacing, corresponding to /3 = 6.0 and /3 — 6.2, with quarks described 

by a mean-field-improved Sheikholeslami-Wohlert action. As a by-product, we obtain the leptonic 
decay constants of B and D mesons. We also present matrix elements relevant for D'' — D'' mixing. 
Our results are summarized in the Introduction. 
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■ The study of Bd—Bd oscillations enables measurement of the magnitude of the poorly known Cabibbo-Kobayashi- 
Maskawa (CKM) matrix element, Vtd, and thus the determination of one of the sides of the unitarity triangle. The 
T-H , frequency of these oscillations is given by the mass difference, 

o • 

O ; Amd = - , (1) 

' ^ ^ where Mg^ and are the heavy and light mass eigenvalues of the Bd—Bd system. Ato, is experimentally mea- 
' surable from tagged Bd meson samples. ^ It is also calculable in the standard model. Keeping only dimension six 
operators after an operator product expansion in which the top quark and W boson are integrated out, the standard 
(— I \ model prediction for Am, is, to next to leading order (NLO) |^,^: 

• ^ ■ A.-n - A//-2 IT/ T/*|2„ c- \r \{^d\Od^='^{^J')\Bd)\ 

Am, - — I VtrfVjf,! r]Bbo(xt)CB(p) , (^j 

d ' where Xt = m^/M^, So{xt) ~ 0.784 Xj-^^ (to better than 1%) is the relevant Inami-Lim function S, the renormal- 
ization scale, 0^^^^ the four-quark operator [57'^(1 — 7^)rf] [^7/^(1 — 7^)*^] and 7]b = 0.55 and Csf/i), short-distance 
coefficients. The renormalization-scale dependence of CsifJ-) and of the hadronic matrix element cancel such that 
Arrid is /^-independent, to the order in perturbation theory at which CsifJ-) is calculated. In the naive dimensional 
regularization modified minimal subtraction (NDR-MS) scheme assumed here. 



J ^^^^ C^) 
^5 - ^ . (3) 



Since |Vtf,| is equal to unity to very good accuracy, a measurement of Am, clearly enables the determination of \Vtd\- 
The accuracy of this determination is limited, at present, by the theoretical uncertainty in the calculation of the 
nonperturbative, strong- interaction effects in the matrix element {Bd\Of^^'^{ii)\Bd). 

An alternative approach, in which many theoretical uncertainties cancel, is to consider the ratio, Ams/ Arrid, where 
Arris is the mass difference in the neutral Bg-Bg system. In the standard model, 



^For a recent experimental review see, for instance, |l| or 
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where 0^b=2 ^j^g gg^j^^g operator as Of^^^, with d replaced by s, and where we have omitted the renormahzation- 
scale dependence, as it cancels in the ratio. Because the unitarity of the CKM matrix implies |Vts| — jVcbl to a few 
percent and a clean extraction of |Vch| can be achieved by analyzing semileptonic B decays a measurement 
of Arris/ Amd yields a determination of \Vtd\- The high frequency of Bs~Bs oscillations makes this a challenging 
measurement. Nevertheless, the experimental lower bounds obtained on Arris [Q,^ already yield interesting constraints 
on the unitarity triangle [p],|7|-|l^. 



The matrix elements that appear in Eq. are traditionally parametrized as 



B,\Of^='{^)\B,) = -M'sJ's^BbM 



(5) 



where q = d or s and where the parameter, Bb^{p), measures deviations from vacuum saturation, corresponding to 
Bb^ (/i) = 1 . Here, Jb^ is the decay constant deiined by 



{0\h^^''q\Bg{p)) = tp^fE 



(6) 



One also usually introduces a renormalization-group invariant and scheme-independent parameter, Bb , which to 
NLO in QCD is given by 

B^': = Cb{p)BbM ^ (7) 

where CBin) is given by Eq. (||) as long as Bb^(p) is computed in the NDR-MS scheme with five active quarks. For 
consistency with the value of rjB given after Eq. (||), should be taken to have its two- loop value with = 225 MeV. 

In the present paper we report on high statistics, quenched lattice QCD calculations of matrix elements and the 
corresponding 5'?7(3)-breaking ratios relevant for neutral _B-meson mixing. We obtain the 5L/(3)-breaking ratio rsd in 
two ways: (1) by calculating Bb^/Bb^ and /b,/ fs^ ^^id combining these two ratios with the experimental mass ratio, 
MbJMb^ ("indirect" method); (2) by calculating the matrix elements, {Bd,s\Off^'^\Bd^s)^ directly and taking their 
ratio ("direct" method), as suggested in Ref. [l3| 0. They mainly differ in the required light- and heavy-quark-mass 
interpolations and extrapolations, since for the "direct" method it is the matrix element and corresponding ratio 
that are interpolated and extrapolated, while for the "indirect" method it is the i?-parameters, decay constants, and 
corresponding ratios. 

As described in more detail in Section |l[ these calculations are performed at two values of the lattice spacing, a 
(~ (2.0 GeV)~^ and ~ (2.7 GeV)""'^), with relativistic Wilson fermions. In order to keep discretization errors in check, 
the lattice calculation is performed with heavy quarks whose masses are around that of the charm and the results 
are extrapolated to the mass of the h. Even in the charm sector, however, quarks have compton wavelengths that 
are not much larger than our lattice spacings, and it is important to reduce discretization errors as much as possible. 
We attempt to do so by describing quarks with mean- field- improved [Q, Sheikholeslami-Wohlert (SW) actions [ p^ . 
When combined with improved operators, these actions lead to discretization errors, which are formally smaller than 
those generated with an unimproved Wilson action [O {aso) instead of O (a)), and which may be numerically smaller 
than those brought about by a tree-level-improved SW action. It is important to note that, as far as four-quark 
operators are concerned, nonperturbative O (a)-improvement has not yet been undertaken, and all lattice calculations 
of Bd^s—Bd^s mixing matrix elements have, as we do, O (a^a) discretization errors, or worse, O (a) errors. It is also 
important to remember that O (a^) errors can be significant in the presence of heavy quarks in a relativistic approach. 

Alternatively, one could take an effective theory approach and work with static, nonrelativistic QCD (NRQCD) or 
Fermilab quarks. In these approaches, the matrix elements are expanded in inverse powers of the heavy-quark mass 
to remove it from the long-distance dynamics. One important advantage is that discretization errors are no longer 



^Our method differs slightly from that of Ref. in that we actually calculate {{Bs\0^'^=^\Bs) / (BalOf^^^lBd)) x 
[Mb^/ Mbs) and multiply by the experimentally measured value of Mb^ / Mb^^- 



2 



enhanced by this mass. An accurate description of the physics of the b quark, however, requires one to consider cor- 
rections in inverse powers of rrib and the calculation of these corrections is made difficult by contributions proportional 
to inverse powers of a. The effective theory and relativistic approaches should thus be viewed as complementary. 

An additional feature of our calculation is that we extrapolate 5t/(3)-breaking ratios in the heavy-quark mass 
directly instead of first extrapolating numerator and denominator and then taking their ratio. The heavy-quark-mass 
dependence will cancel partially between numerator and denominator and therefore make the extrapolation more 
reliable. This approach turns out to be particularly fruitful for the determination of fs^/fB, where the statistical 
error is significantly reduced by a direct extrapolation of the ratio. 



Our main results are 



f^^V^ = 1.15(2) , r..= (^'^'^'-^ 



fB,^: - 210(21) til MeV , g^V^ = 0.89(7) 
= 241(14) tl' MeV , g;^^ = ^'^^^^^ 
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1.38(3) 



+0 



BbJBb,= 0.98(2) II 



where the first error is statistical and the second corresponds to the systematic uncertainties added in quadrature. In 
quantities involving ratios of B parameters, the renormalization-scale dependence is not specified, as it cancels. We 
consider a wide array of systematic uncertainties, as discussed in Section VII. Wc normalize dimcnsionful quantities, 
involving decay constants, by /u^ because some systematic (and statistical) uncertainties, including possibly those 
associated with quenching, partially cancel in the ratio. The original quantities can then be recovered by using the 
experimental measurement of fu. 



The results of Eq. (jg) will be compared, in Section |VII|, to earlier calculations of some or all of these quantities 
performed with propagating heavy quarks |l^,|2^,|2^ and with non-relativistic quarks [||,|26). The comparison of the 
B parameters with results obtained using static heavy quarks p7|-|30|] will be addressed elsewhere |^ . 

Because decay constants and the corresponding matrix elements are necessary for obtaining the results of Eq. (H), 
we also have results for these decay constants. We find 
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nun) tf^ Mcv 
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fu - 210(10) MeV 
fo, = 236(8) tZ McV , 
^ - 1-13(2) tl , 



(9) 



where the first error is statistical and the seco nd sy stematic, as discussed in Section VII. A comparison with recent 
quenched results |32 -3^ will be made in Section VIII and a discussion of unqucnched results |^-El| v 



will be undertaken 



^The results of preliminary analyses on the same lattices were presented in Refs. |]l^-|l8t. 
*We take Mb, = 5279 MeV, Mb, = 5375 MeV, Md = 1864 MeV and Afo, = 1969 MeV. 



^Ref. gives fn = 241(32) MeV as a summary number, which is in good agreement with our determination of this 
quantity (see Eq. (BI) below). However, recent determinations appear to yield larger values, albeit with large uncertainties: 



fn, = 285(20)(40) MeV |||, fn, = 323(44)(36) MeV H and fo, = 280(19)(44) MeV [||, 
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in Section VII F . Note that results for fs can be combined with the measurement of the branching ratio for the rare 
decay t^>^t, when it becomes available, to yield a clean determination of \Vub\ - § 

While short-distance D'^—D^ mixing is highly suppressed in the standard model [Q, it can be enhanced in super- 
symmetric extensions p9|-|53||, above even the long-distance contributions discussed in Ref. |Q. Thus, we give the B 
parameter and decay-constant combinations relevant for the matrix element of the left-left, AC = 2 operator, which 
is one of the operators that can contribute in supersymmetric extensions: 



Si5(Md) =0.82(3) 



222(10) 
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where is obtained by multiplying Bd{Md) by CoiMo), with 



An 



Ji = 



6719 



3750' 



(11) 



where the two-loop as{fJ.) is evaluated with Ai-i = 350 MeV 



II. SIMULATION DETAILS 



Our results are based on quenched, SU{3) gauge configurations, calculated on a 24'^ x 48 lattice at /3 = 6.2 and a 
16^^ X 48 lattice at /3 — 6.0. The configurations are generated using the hybrid over-relaxed algorithm described in 
Ref. [ p5[ . The parameters of the simulations are summarized in Table |. 

We describe quarks with the Sheikholeslami-Wohlert (SW) action |^ ^ 



Csw " 



(12) 



where is the standard Wilson action, go the bare gauge coupling, P^^, a lattice definition of the field strength 
tensor, the appropriate quark hopping parameter and Cswi the so-called clover coefficient. Here, ip stands for both 
light (q) and heavy (Q) quarks. While with the Wilson action (csw = 0) spectral quantities suffer from discretization 
errors of 0{a), the tree-level value, Csw — 1, guarantees that these errors are reduced to 0{asa) ||l^,Q. In the 
present paper we work with a mean-field estimate of the clover coefficient [^4[Q, = 1/mq with uq — (■|TrJ7p)3. 
Since this estimate accounts for large tadpole contributions and is closer to the nonperturbative value of the clover 
coefficient | ]58| which removes 0{a) errors to all orders in ag, our discretization errors may be numerically smaller 
than for Cgw = 1- It is important to note, however, that in the presence of heavy quarks with masses amg ^ 0.3 or 
more, discretization errors of O {^{arnqY) may be comparable to those of O (asOTOg). 

In order to improve matrix elements up to O (aso), we must further "rotate" the quark fields [^6|| : 



1 - 



Z7 • - (1 - z)m^ ))''/') 



(13) 



where is the symmetric covariant derivative, the bare quark mass to be defined later and z is a real parameter 
which can have any value between and 1. For heavy quarks, a large source of discretization errors is the mismatch 
in the normalization of tree-level, zero-momentum, continuum and lattice quark propagators. Eq. (^3|) with 2 = 
corrects this mismatch at O (a) . Since we can compensate this mismatch completely by implementing the El-Khadra- 
Kronfeld-Mackenzie (EKM) normalization p9||. 



we choose the latter instead of Eq. m 



(14) 



For recent reviews of lattice calculations of the quantities presented in Eqs. (g) and (|9|), please see 
''Here and below tj^^ — | [7^1 Ti-], where 7^ are the usual set of Euclidean 7-matrices with {7,^,71^} = 2(5^ 
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At both values of the lattice spacing, we work with several values of the heavy-quark hopping parameter straddling 
the charm. This enables us to extrapolate our results in heavy-quark mass from the charm sector, where discretization 
errors appear to be only a fraction of the result, to the bottom sector, where these errors would be very large were 
we to perform the simulation directly with such quarks. We also consider several values of the light-quark hopping 
parameter around that of the strange. Then, we interpolate our results in the light-quark mass to the strange and 
extrapolate them to the chiral limit. The values of the hopping parameters used in our paper are given in Table ||. 
For completeness, the masses of the corresponding light-light and heavy-light pseudoscalar mesons are also given in 
physical units. 

To isolate the ground state more efficiently in the correlation functions that we calculate, we use fuzzed sources 
and/or sinks |]60| . These are extended interpolation operators that have improved overlap with the ground state. Of 
course, operators whose matrix elements we wish to compute are kept local. 

Statistical errors are estimated from a bootstrap procedure which involves the creation of 1000 bootstrap 
samples from our set of 188 (498) configurations at (3 — 6.2 {(3 = 6.0). Correlators are fitted for each sample by 
minimizing x^. The quoted statistical errors are obtained from the central 68% of the corresponding bootstrap 
distribution. 

To convert our values for decay constants into physical units we need an estimate of the inverse lattice spacing. 
The determination of this quantity is discussed in Section 



III. MATCHING AND RUNNING 



Because the lattice and the continuum treat ultraviolet modes differently, the extraction of continuum matrix ele- 
ments from lattice calculations requires a matching procedure. Ideally, this matching is performed nonperturbatively. 
For our mean-field-improved action, however, nonperturbative matching coefficients are not available and we resort 
to perturbation theory instead. 

The simulation is performed with the fermion action of Eq. ( |l2|) with Csw = 1/^*0 standard Wilson gauge 

action. However, the fully mean-field improved action would involve normalizing each occurrence of a link variable 
in the action by the measured value of uq. To recover the results we would have obtained had we used the latter, 
we must, in interpreting the results of our simulation, use rescaled bare couplings cXg = Os/uq and = k^uq, with 
Us = (7Q/(47r). already has its desired rescaled value and no additional rescaling is necessary, is actually only 
a first guess at an improved expansion parameter and one may try to optimize this choice [ p^ . This issue will be 
elaborated on at the end of the present section and, for the moment, we will generically denote the coupling by a^. 

In perturbation theory, the effect of normalizing link variables is obtained by expanding the factors of uq in powers 
of the strong coupling. At O (as) we have jl^ , 

/I \ ^ a 47r^ 

uo = (-TTUpj =1 + ^X withX = -— . (15) 

This means that every occurrence of in a first order perturbative expression must be replaced by {l~{as/4:Tr)X). 
Because factors of Cgw always appear multiplied by at least one power of in perturbative expressions, Cgw = 
1 — 3(as/47r)X can be replaced by Cgw — 1 hi first order expressions. This is what we do to determine the central 
values of our results. However, in obtaining errors, we consider the variation induced by taking Csw = l/'^o- 

To extract the B parameters of Eq. from ratios of three-point and two-point correlation functions on the lattice, 
we must match the mean-field improved, EKM-normalized lattice axial-vector current to its Euclidean continuum 
counterpart via 

^ ZA{as)A'^\a) , (16) 

with 

Ajf*(a) = + arhg ^1 + arfiQ ^2Rq 2Rq (Qj^-fsq) (a) , (17) 
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where arUq^Q = {^/nq^Q — l/Kcr)/2, with the mean-field improved version of the critical hopping parameter, Kcr, 



which is determined nonperturbatively as detailed in Section 
loop|: 



Using the results of Refs. |63, p3|, we find at one 



ZA{as) = 1 



47r 



7m 75 



As,) -X 



1 - -^[7.90 + 0.33csv 
An 



3.00c^„] , 



(18) 



where, in the notation of Ref. |62[ , A^^-^^ and A^^ arise from the one-loop corrections to the vertex 7,^75 and to the 
quark wave- function, respectively. The effect of mean-field improvement is encoded in the term proportional to X . 
Without mean-field improvement, i.e. X = 0, the coefficient of Q!s/47r would be substantially larger: 18.39 instead of 
5.23 for Csw = 1- 

The matching of the four-quark operator O^^^^ {F stands for the flavor of the heavy quark) is complicated by 
the fact that Wilson-type fermions break chiral symmetry explicitly, inducing mixing amongst four-quark operators 
of different chirality. The following five operators form a complete basis for this mixing on the lattice in the parity- 
conserving sector: 



= 7m X 7m ± 7/^75 X 7^75 
O^^* = / X / ± 75 X 75 , 



(19) 



where a sum over Lorentz indices is implicit and where F x F stands for (1/a^) (l-}-am„) {1 + amn) (2K„2Kn) 



{QTq){QTq). The parity even component of 0^^~^ corresponds to the continuum equivalent of so that, at one 



loop: 



Of^=\^i) Zn(a„ a/i) Ol'^* (a) + ^ZH(a,)Ol'^' (a) 



(20) 



i=2 



Zii has a logarithmic dependence on aji, where jj, is the scale at which the continuum operator is rcnormalizcd, while 
Zij, i = 2j_---,5, remains finite as a vanishes. Mean-field improving the dimensional reduction (DRED) results of 
Refs. fs^-p^ and using the matching between DRED and NDR-MS given in Ref. |^, we find that one-loop matching 
to the NDR-MS scheme is given by: 



Zii{as,afi) = 1 + 

47r 



41n(aAi) + i{5(A^^ 



*7m75 



-(A/ + A^J + 8Asi} - 2-2X 



1 ^[41n(a^) + 24.52 - 9.33cs 

An 



(21) 



4.88cf 



■^12 (as) 
■^13 (as) 



11a,.. . - a, - 
--^[A.-A,J=-^[8.84-9.15cs 



3.13d 



4 2 2 

- — Zi2{(^s) , Zii{as) = —Zi2{as) , ^i5(as) = — ^i2("s) ' 



where the Ar, F = 7^, 7^75, / and 75 arise from the one- loop corrections to the bilinear vertices associated with F. ^ 

For the numerical evaluation of the renormalization constants, we choose to work with the MS coupling, a^, 
obtained via 



^5^(3.41/a) = av(e^/S3.41/a) (l + 2^) , 



(22) 



^In practice, we use where results are given for arbitrary Csw and where loop integrals are calculated to higher numerical 
accuracy. 

^Here again, we use the more precise and general Csw results of Ref. [M for Ar. 
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solving 



where ay is the coupling defined from the heavy-quark potential. The latter is obtained from our simulations by 

In(uo) = ^ay(3.41/a) (1 - 1.185ay) . (23) 



Values of a-^ for scales different from 3.41/a are obtained by solving the two-loop running equation numerically, with 
Uf — Q. Both a-^ and ay have been shown to lead to expansions that are much more convergent than those in terms 



of the bare lattice coupling [g^. For completeness, we give values of a-^{ix) for a collection of /i in Table [II 



Having chosen the coupling, we must fix the scale, q* , at which it is evaluated. We take 2/a as a central value but 
allow q* to vary from 1/a to it /a to estimate the uncertainty associated with this choice. This range for q* covers 
typical, ultraviolet lattice scales. 

As it is convenient for the heavy-quark extrapolations and does not generate large logarithms, we match the lattice 
results for the B parameters at = Mp, , where Afp, is the mass of the heavy-light meson composed of a heavy 



antiquark Q and a massless quark, I — d or u. As we will see in Section Vl, the heavy-quark extrapolation yields 
{^Bd ) ■ Values for the renormalization-group- invariant and scheme-independent B parameters are then simply 
obtained from Eqs. (|) and (|): B'^° = CBiMB^BB^iMBj- Values for with ^ ^ Mp^, can also be obtained 

straightforwardly through: i?p^(/i) = [Cp(Afp^)/Cp(/i)] i3p^(A/p^), with or without an expansion in as of the term 
proportional to J5 in Eq. (|^) 

Since we match the matrix elements and B parameters defined in Eq. (|^) at one loop, we may choose to expand 
combinations of renormalization constants such as Z\\IZ\ or ZnZu, j = 2, • • • , 5, to order a^. In the present paper, 
all central values are obtained without expansion. Nevertheless, we have checked that expanding these combinations 
makes negligible differences in the final results. This is because the one-loop corrections to the renormalization 
constants are small, especially after mean- field improvement. 

IV. CORRELATION FUNCTIONS 

To determine IS.F — 2 matrix elements and their B parameters, as well as decay constants, we compute the following 
two- and three-point functions: 

C'A'p{t.) = Y,{Pl{x)A'S\Q)) , (24) 

X 

S,y 

as well as Cp^{tx) and Cpp{tx)- Here, the subscripts F on operators indicate that they are fuzzed |^^, while operators 
with no F are local. For the correlation functions, the superscripts F and L indicate which of their operators are 
fuzzed (F) or local (L). In Eq. (p^), A{f* is the time component of the EKM- normalized axial- vector current defined 



in Eq. (17), Pp is a fuzzed [|60[ version of the pseudoscalar density Qj5q, and O stands for any of the four-fcrmion 
operators in Eqs. ( p^ and (|20|). In the present section, we set the lattice spacing a — 1 and omit light-quark indices 
on heavy-light quantities, for notational simplicity. 

At large Euclidean time separations, the three-point correlation function Cq'"^ {tx,ty) has the asymptotic behav- 
iorQ 

C%^P(t,,ty) ^ ^{0\Pl\P){P\O\P){P\Pl\0)e^p[-Mp{T -t, +ty)] , (25) 



^"On our periodic lattices this corresponds to aty and a{T — tx) much greater than 1 but small enough so that the desired 
time ordering dominates. Here, aT = 48 is the time extent of our lattices. We consider times and ty such that T/2 <tx<T 
and Q <ty < T/2. 
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where Mp is the mass of the pseudoscalar P. Therefore, fits to the ratios of correlation functions: 



fiPLF (f f \ 

C\p{tx)C\p{ty) 



i;ZiBp , (26) 



^ ^ ' (27) 



C'p^{tx)C^'p{ty) {{)\Pf\P){P\Pf\0) ' 



where Bp is the B parameter corresponding to the heavy-hght meson P, yield the desired quantities up to renormal- 
ization constants that we determine perturbatively (see SectionElS) and factors of Mp and of the type, (0|Pf|P 



that we determine from fits to the two-point functions CpJ;(ia:)|^ and Cpp{tx), respectively. Note that in Eqs. ( [25| ) 
and (p7| ), all pseudoscalar meson states have vanishing three momentum. An example of a plateau for the ratio of 
Eq. ( p6| ) is shown in Figure 

To determine the decay constants, we consider 



C'p^tx) ZA{0\Pm 



tanh[Mp(T/2-t^)] . (28) 



In order to investigate the dependence of the matrix elements and B parameters on heavy- and light-quark mass, 
we need the heavy-light and the light-light pseudoscalar meson masses, and the light-light decay constants, which we 
use to set the scale. As already stated, we obtain the heavy-light pseudoscalar meson mass, Mp, from a fit to the 
large time behavior of Cpp{tx). Similarly, we get the light-light pseudoscalar meson mass, mp, from a study of the 
two-point function CpJ"{tx) where p — Qi75(j2 is the pseudoscalar bilinear made from the two light quarks qi and g2- 
Finally, the light- light pseudoscalar decay constant is obtained from a fit to the ratio CpJ^{tx)/Cpjf{tx), where ag is 
the EKM- normalized version of 917075Q2, much in the same way the heavy-light decay constant is obtained. 



The time ranges over which the various correlation functions and ratios are fitted are given in Table IV. 



V. LIGHT-QUARK-MASS EXTRAPOLATIONS AND DETERMINATION OF THE LATTICE SPACING 

Results for physical Bd and Bs mesons require investigations of the dependence of the lattice measurements on 
light- and heavy-quark masses. We begin by the light-quark-mass extrapolations and interpolations. To obtain K^r, 
the critical value of the quark hopping parameter, we study the behavior of the light-light pseudoscalar meson mass, 
rup, as a function of ki and the hopping parameters of the light quarks which compose it. We assume that it 
obeys the partially conserved axial-current (PCAC) relation 

(oTOp)^ (ki, K2) = Pmiarhi + 07712) , (29) 

where is the 0(a)-improved quark mass, given by 

rhi = mi{l + bmarrii) , (30) 

with ami = 1/ (2^.;) — l/(2Kcr) | ]6^ . At tree level, which is sufficient with our mean-field-improved action, — —1/2. 
Kcr is then obtained by fitting the function of Eq. ( p9|) to our results for (am^)^, with and Kcr as fit parameters. 
Ki and K2 are taken amongst the values in the row labeled Kq in Table |[ Addition of a term quadratic in quark mass 
in Eq. (^^ makes very little difference to the central value for Kcr and this variation has a negligible effect on the 
light-quark-mass extrapolations of the matrix elements of interest. The linear extrapolation of Eq. (|2^) is shown in 
Figure ||. 

For mp and all the quantities we study here, in addition to higher-order polynomial corrections in light-quark 
mass, chiral perturbation theory also predicts the presence of chiral logarithms. These logarithms, however, are 
difficult to isolate numerically and are modified by the quenched approximation [p7H69|. Thus, we perform polynomial 



^Cpp gives a particularly good signal for Mp. 
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interpolations and extrapolations from our intermediate values of light-quark mass, only considering in most cases 
the leading, non-trivial dependence on light-quark mass. 

Next we determine the lattice scale and the bare strange-quark mass, arhs, in lattice units. We obtain both 
quantities simultaneously from the kaon's decay constant, /k, and its mass, rriK- As long as the mass dependence of 
fp, the decay constant of a light-light pseudoscalar meson, can be described by a function of rrip only, Eq. ( pof ) implies 
that both quantities depend only on the sum of quark masses and that only an interpolation in these masses, not an 
extrapolation, is needed to obtain afx- The use of /k, instead of, for instance, rriK' , is prompted by the fact that the 
only quantities we report on here which depend strongly (i.e. not logarithmically) on are the heavy-light decay 
constants. Our procedure has the added benefit that it is also applicable in unquenched simulations |^ Thus, we fit 
our results for aj^jZ a to: 

= a/ +/3/(a'7ip)2 -|-7/(aTOp)'^ H . (31) 

We find that this parametrization describes our results well. Assuming that the lattice spacing is fixed with Jk (or 
equivalently mpc), we solve Eq. ( ^l[ ) for afp (or arrip), at the point specified by the physical ratio mpc/fK, with 
fx = 159.8 MeV and niK = 493.7 MeV. arhs is then obtained from the resulting amji, using our earlier fit to 
Eq. (p9|), and the lattice spacing from the resulting uJk (or arriK)- The fits at our two values of the couphng are 
shown in Figure H. Because of the slight curvature, we favor the quadratic fits. These fits also give a value of /k/Ztt 
which is closer to the experimental result of 1.22: 1.16(3) (1.19(2)) for quadratic fits instead of 1.131(10) (1.147(7)) 
for linear fits, at /? = 6.2 (6.0). In any case, linear and quadratic fits give nearly identical results for the scale and 
aifis- The values of Za used are those obtained from Eq. ( [l^ ) with Ug = aj^{2/a). Systematic uncertainties in the 
dete rmination of a^^ and afhs will be addressed when we discuss the uncertainties on our B-physics results in Section 
VII . A summary of the results for a~^, Kcr and afhs used below is given in Table 0. 



Having determined the strange quark mass and critical hopping parameter, we interpolate and extrapolate our 
various heavy-light matrix elements in light-quark mass to these values. We assume that the up and the down quarks 
are massless. This is an excellent approximation for the quantities we study. As mentioned above, we perform 
polynomial interpolations and extrapolations from our values of light-quark mass. Thus, we fit all quantities of 
interest, Y , in lattice units, to the functional form: 

Y{KQ,Hq) ^ ay + |3Yam,q^ . (32) 

In Figures |5[ ^ and we exhibit the light-quark-mass dependence and the corresponding fits to Eq. ( ^^ for the 
quantities l/{aMp^), Bp^ , {Pq\0^^^'^\Pq) and afp^/ZA^ for all values of heavy-quark mass at /3 = 6.0 and 6.2. For 
clarity of presentation, the B parameters and {Pq\0'^^^'^\Pq) arc renormalized at a common scale of 5 GeV |^. In 
all cases, the light-quark-mass dependence is mild and, to good accuracy, linear. However, we do observe that the 
matrix elements of O^^^^ have a stronger dependence on light-quark mass than the other quantities. 



VI. HEAVY-QUARK-MASS EXTRAPOLATIONS 



The second extrapolation we have to perform is in heavy-quark mass. We use heavy quark effective theory (HQET) 
as a guide, with Mp^ as a measure of the heavy-quark mass everywhere except for kinematical dependencies, where 
the appropriate meson mass is used. As before, / stands for either a d or a u quark. Other choices for the heavy-quark 
scale, such as Mp^ , for instance, make very little difference to the final results. We match our QCD results onto HQET 
at the heavy-quark scale A/p, , and cancel the leading logarithmic dependence on Mp by including terms of the form 
[as(Afp, )]~^o /2/3o^ where 70*^^^ are the relevant, one-loop anomalous dimensions in HQET, and /3o is the one-loop 
/3-function coefhcient. 7o'^'^^ = —4 and —8 for the decay constant and four-quark matrix element, respectively ||70|,[7l|. 
Thus, we define 



^■^The K* is not a stable particle once light-quark loops are allowed. 

'^^These values of Bp^ and a'''{Pq\0^^^^\Pq) are obtained by matching the lattice results onto the NDR-MS scheme at the 
scale 2/a, then running to 5 GeV in the MS scheme at the two- loop level with Uf — and the coupling constant described in 



Section 



III. 
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for the decay constants, 



for the B parameters, and 



Za 



2//3o 



$|(Mp,)^i?P,(Afp,)xK(Afp,)]«/'3o , 



(33) 



(34) 



'i'ip=2 



AP=2 



aM 



(35) 



for the AF — 2 matrix element. Similar scaling functions are defined for S'C/(3)-breaking ratios, in which the leading 
logarithmic dependence on heavy-quark mass cancels. In Eqs. (^), (Q) and (|35|), we evaluate as{M) through 
27r/(/3o log(Af/AQCD)) with Aqcd = 100 MeV and (3o = 11 — (2/3)n/ with n/ = 0, since we are working in the 
quenched approximation. This one-loop coupling approximates the lattice couplings defined through Eqs. (|2^) and 
(2^) rather well, for the values of M required here. In fact, final results depend weakly on the value Aqcd used in 
the heavy-quark extrapolations (see also Figure || and discussion below) . 

For X{Mp^) = $AF=2 (-^^fl ) ' '^'b(-^^Pi)' ^'fi^^Pi) s-iid the corresponding scaling functions for S'C/(3)-breaking ratios, 
we use the HQET-inspired relation 

X{Mp) = Ax 1 + Bx(^^]+Cx 



aMf 



(36) 



to investigate the heavy-quark-mass scaling of these quantities. The leading logarithms make little difference in the 
extrapolation, as shown in Figure ^ where we plot the extrapolation of the AF — 2 matrix element, which has the 
strongest logarithmic dependence amongst the quantities we study, with and without these logarithms. 

In Figures ^ and |l^ we display our results for $^(Mp) and *I'af=2(-^^J^; )' constructed from Bp^{Mp^) and 
{Pq\0^^^'^ {Mp^)\Pq) renormalized at Mp in the NDR-MS scheme, as functions of inverse heavy-meson mass at 



/? = 6.0 and 6.2 and for q = s and I. In Figure 0, we plot <^){Mp,) versus l/(aMp, ). Finally, in Figures |l|,|l3|and|l^ 
we plot the extrapolations of the corresponding SU (3)-breaking ratios. The fit parameters of the heavy-quark-mass 
scaling of the various quantities studied are summarized in Table VI. While the heavy-quark-mass dependence of 



most quantities is mild, that of $y-(AfpJ and especially ^'AF=2(-^^fi) quite severe. 

We extrapolate SU (3)-brcaking ratios directly in heavy-quark mass because the mass dependence cancels partially 
between numerator and denominator, making the extrapolation less pronounced and, thus, more reliable. This is 
especially visible for fp^ / /p, (Figure ^4|) , where the extrapolation of the ratio yields much smaller uncertainties than 
the ratio of the extrapolations. In all cases, the extrapolation of the ratio is in excellent agreement with the ratio of 
extrapolations. It is interesting to note, also, that the heavy-quark-mass dependence of /p, //p, appears to be the 
same as that of MpjMp^. 



VII. SYSTEMATIC UNCERTAINTIES 



and [X 



In these tables, 

the first error on each quantity is statistical. The remaining uncertainties are systematic and we discuss them now. 



Our main results at the two values of lattice spacing are summarized in Tables VII , VIII 



A. Discretization errors 



Ideally, one would extrapolate all computed quantities to the continuum limit, where discretization errors vanish. 
With two values of the lattice spacing, however, this is not possible. We must therefore use the information that we 
have to estimate the uncertainty associated with residual discretization effects. 



In Table VII 



results for the decay constants display some dependence on lattice spacing. This suggests that 
discretization errors for these quantities may be important. The leading discretization errors with the mean- field- 
improved SW action are formally of O (a^a), as they are for the tree-level improved SW action. Subleading errors 
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begin at O (a^a) |^ and O (a^) . To estimate these leading and subleading errors, we consider the following variations 
in our procedure. O (a^a) improvement of the axial current requires one to include the effect of the adf^P counterterm 
through the replacement (P is the pseudoscalar density) 



as well as to rescale the quark fields as 



-^A^ + CAad^P , (37) 



i/) ^ (1 + ^am^)^p , (38) 



with both ca and bA evaluated at one loop [|66|,[72[. From a comparison of results obtained with ca and bA set to their 
tree-level values {ca — and bA ~ 1) to those obtained with ca and bA evaluated at one loop, we can estimate the 
effect of O (agfl) discretization errors. We do not use the one-loop results as central values for the decay constants to 
be consistent with our determination of the B parameters. Indeed, O (ago) improvement of the four-quark operators 
would require one to consider the mixing of these operators with operators of dimension seven, which is beyond the 
scope of the present paper. 

As already mentioned, to subtract higher-order discretization effects, we use the EKM normalization of Eq. (|lj). 
Thus we define the one-loop variation in the normalization of the quark fields through 

where, using the results of fz^ , we find b]^^°°^ = 1 -f (as/47r) [24.03 + X] = 1 + 10.87(as/47r) for Csw = 1- c]^^°°^ 
is given by — 1.20(as/47r) | ]73| . We find that the replacements of Eqs. ( [37| ) and (|3|) have opposite effects. The 
former lowers the decay constants while the latter increases their values. For the one- loop values of ca and bA, the 
cancellation is rather good and the resulting one-loop versus tree-level variation is certainly an underestimate of the 
residual discretization errors proportional to a. To get a more realistic estimate, we consider the variation brought 
about by each replacement separately. These variations are shown as the second error on the decay constants and their 
S'J7(3)-breaking ratios, foj Id and /bJ/b- We take the largest of the two to be a measure of residual discretization 
errors proportional to a. Fj 

To get a handle on errors proportional to higher powers of a, we consider the result of using the tree-level, quark- 
field normalization of Eq. ( |3^ ) with bA = I, instead of the EKM normalization of Eq. (p^). These two normalizations 
differ at O (a^) and we take the resulting variation to be a measure of these additional discretization errors. This 
variation is shown as the third error on the decay constants and their S'[/(3)-breaking ratios, fos/fo and JbsI fs- 

The symmetric discretization error that enters the systematic error in the final results of Eq. (|^) is obtained 
by combining in quadrature our estimates of the residual discretization uncertainties proportional to a and of the 
uncertainties proportional to higher powers of a. While these two uncertainties are comparable at /3 = 6.2, the latter 
are significantly larger at /? = 6.0 in the 6-quark sector. 

A similar estimate of discretization errors can be carried out for the four-quark matrix elements and their B 
parameters. However, as we have already mentioned a full quantification of O (a^a) effects for these quantities is 
beyond the scope of this paper. In fact, many discretization effects, such as those associated with the normalization 
of quark fields, cancel or partially cancel in the ratios of matrix elements used to define the B parameters and SU (3)- 



breaking ratios. Furthermore, in Table VIII, results for B parameters and S'C/(3)-breaking ratios exhibit very little 
lattice-spacing dependence, supporting the idea that discretization errors for these quantities are small. Thus, we 
assume that their statistical uncertainties encompass possible residual discretization errors. For the quantities in 



Table [X, however, which are obtained using the decay constants, we take into account the discretization errors on 



these constants. 



^''O (osa) errors, as well as all errors proportional to a, are absent in non-perturbatively O (a)-improved calculations. 
'^^The largest variation is the one of Eq. (pj]). 
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B. Matching uncertainties 



As already indicated in Section III, to estimate the systematic errors arising from the pertm'bative matching of 
the various quantities we compute, we vary the scale, q*, at which a-^ is evaluated, in the range between 1/a and 
7r/a, and compare this with the result obtained for q* = 2/ a. We also consider the variation coming from computing 
Zii (Mp,), Zii and Za with the constant Csw set to its mean- field- improved value instead of 1, keeping q* — 2/a fixed. 
5t/(3)-breaking ratios of decay constants are not affected by these variations while those of the B parameters are not 



significantly so. These variations are reflected in the fourth error in Tables VII, VIII and [X. 



C. Heavy-quark-mass extrapolations 

As shown in Figure |l|, the decay constants have a pronounced extrapolation in heavy-quark mass, and the term 
quadratic in 1/ (aMp^ ) on the right-habd side of Eq. ( p6| ) contributes significantly. To quantify the systematic error 
associated with this extrapolation-the fifth error on the decay constants-we perform a fit of the heaviest three points 



in Figure 11 to the right-hand side of Eq. (pw, without the quadratic term. For the S'?7(3)-breaking ratios of decay 



constants, we perform a constant fit to these same three points. These uncertainties are propagated to the results of 



Table IX. 



Figure ^ indicates that the heavy-quark-mass dependence of the B parameters and their SU (3) breaking ratio is 
mild and to very good approximation linear. We have verified that a linear fit to the three heaviest points gives results 
that are well within the errors bars of the fit to all five points at /3 = 6.0. We assume that the same would be true at 
(3 = 6.2 if we also had five heavy quarks, as there is no evidence for curvature on the three points that we have. 

The AF — 2 matrix elements have a very pronounced dependence on heavy-quark mass, as seen in Figure |l^. Since 
we are not reporting results for the four-quark matrix elements themselves, we do not quantify the systematic errors 
associated with their determination. 

One may worry that we have only three heavy quarks at /3 = 6.2 in our calculation of AF — 2 matrix elements. 
However, as we have just seen, the heavy-quark scaling of the B parameters and their SU (3)-breaking ratio is mild 
and displays no evidence for curvature. This is confirmed by the behavior at (3 — 6.0 where we have five heavy quarks. 
Thus, we believe that our results for B parameters and the derived quantities of Table |^ at f3 — 6.2 are reliable. The 
situation is certainly not as favorable for the AF — 2 matrix elements themselves. For those, there is evidence for 
curvature and our three points at /3 = 6.2 can be thought to yield only a rough estimate. Thus, we do not attempt 
to give a final result for these matrix elements from this procedure. In fact, this strong mass dependence is one of 
the problems that makes a reliable determination of rsd, from the ratio of individually calculated {Bq\0^-^^^\Bq) , 
q = s,d, difficult. The extrapolation of itself, on the other hand, is much milder and the curvature is much 
reduced. Thus, we extrapolate it linearly and verify, at (3 — 6.0 where we have enough points, that the result of a 
quadratic fit, r^^J^^^ = 1.50(17), is entirely compatible. |^ In any event, the final value of Vsd that we quote is that 
given by the "indirect" method, where none of this is a problem. 

Another concern may be that our lightest heavy quark is too light to be in the heavy quark scaling regime. However, 
in the extrapolations used to obtain the results of Eqs. (||) and (^, the points corresponding to this quark are consistent 
with the smooth curves determined by the other points. Furthermore, our heaviest quarks are as massive as those in 
other relativistic calculations (see for instance Ref. j24|). Thus, we are not distorting the heavy-quark extrapolations 
by including these lighter points, nor are we missing information on the heavier-quark end. Finally, as described 
above, we include in our errors, the variation obtained by ignoring our lightest two points, where appropriate. 

Ideally, one would have continuum extrapolations of results such as ours and of the same quantities computed in 
the static limit (corresponding to an infinite-mass heavy quark). Results for the h would then be obtained by an 
interpolation in heavy-quark mass instead of by extrapolation, as they are here. We leave such studies for the future. 



^^To estimate the heavy-quark extrapolation error at /3 = 6.0, we consider the variation due to the removal of the lightest two 
points from the linear extrapolation of Fig. hS. 
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D. Uncertainties in the determination of the lattice spacing 



In quenched calculations, the value of the lattice spacing varies significantly with the quantity used to set the scale. 
This variation may be due, in part, to quenching effects, as well as any other systematic uncertainty that affects the 
quantity used to set the scale. In this paper, we determine the lattice spacing from Jk and rriK, as described in 
Section 0. We then vary the inverse lattice spacing, a^^, by ±7%. This range covers the variations observed in the 
determination of the scale from gluonic or light-hadron spectral quantities, and with the same action and parameters 
as we use Jz^ p^ , as well as the variation due to the uncertainty in the perturbative determination of Za- 

Uncertainties in the lattice spacing will obviously affect the determination of the decay constants and the dimen- 
sionful quantities derived from them. They will also slightly change the length of the heavy-quark-mass extrapolations 
(Figures |l^, |p and [lo| ). Furthermore, they induce a variation of order ±15% in arhs, which we obtain from the mass 
of the kaon, and therefore affect all quantities that depend on this mass. 

In practice, we find that the variation of the lattice spacing discussed above does not induce a significant change 
in the B parameters. However, it does affect all the decay constants and the corresponding SU (3)-breaking ratios as 
well as the quantities in Table IX which a re obtaine d from these constants. All of these observations are reflected in 
the sixth error on the quantities in Tables VII, VIII and [X. 



E. Determination of the light-quark masses 



Another source of systematic error is the uncertainty in the determination of the light-quark masses. In the previous 
section we saw that changing the lattice spacing, while keeping fixed the physical quantity used to set the strange- 
quark mass, induced approximatively a ±15% variation in arhg. One can also imagine doing the reverse, i.e. holding 
the lattice spacing fixed while varying the observable used to set the strange-quark mass. For instance, we could have 
used the K* or meson masses to fix arhg, instead of mpc- Let us denote the resulting values of the strange quark 
mass by arhs{mK*)^ etc. Due to quenching and other systematic effects, the values obtained may differ. 

To estimate these differences, we turn to the literature. In Ref. where the determination of quark masses with 
different fermionic actions is reviewed thoroughly, Bhattacharya and Gupta conclude that for mean-field-improved SW 
actions such as ours, afhsirrifj^) / aihs(rnK) — 1.2, a statement that they find depends very little on the gauge coupling 
/3, at least in a range that covers our two simulations. This result is obtained with the lattice spacing fixed by the mass 
of the p meson. With the lattice spacing set in this way, it is straightforward to show that arhs{rnK*) — arhs{m^). 
This follows from the observed linear behavior of the light vector meson masses and the fact that rup ~ 2m/f* — 
in nature. The picture changes when the value of the lattice spacing differs from the one given by the mass of the 
p, a(mp). In the present case, however, the values of the lattice spacing that we use appear to be consistent, within 
errors, with the values of a{mp) that can be inferred from Ref. ||75| . Thus, we consider that the 20% upward variation 
in arhs described in Ref. p5[ is a reasonable estimate of the uncertainty associated with the different possible choices 
of an observable to fix this mass. The changes that this variation induces in our results are reflected in the seventh 



error on the quantities in Tables VII, VIII and IX 



F. Quenching errors 

Quenching effects for the quantities of interest here have been studied using quenched chiral perturbation theory [ |69[ . 
They are typically a few percent for the B parameters if the theory's couplings are constrained by large- A'c arguments 
and by the reasonable range of Ref. and larger outside these ranges. Recent results for the decay constants, 
obtained with two flavors of dynamical quarks (n/ = 2), show little variation in fB^I Ib compared to its quenched 
value p8|-p[|. The authors of [Q flnd that this ratio is enhanced by (5 ±3)% in their calculation with light dynamical 
quarks. Thus, if quenching effects on the B parameters are small, commensurate variations on r^d and ^ are expected. 
For the decay constants of B (Bg) mesons, quenching effects appear to be significant p8|-^. For instance, the authors 
of Ref. ||3^ find that these decay constants are enhanced by 11% (14%) when light-quark loops are included, with a 



^'^Some of the baryons considered in Ref. would give lattice spacings outside our ±7%. However, these particles are more 
susceptible, than the particles we are studying here, to systematic effects such as those associated with the finite volume of the 
lattice. 
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statistical significance of 2 to 3 standard deviations. For the D (Ds) mesons, the effect is of 3% (7%) and consistent 
with zero within roughly one standard deviation. These latter results suggest that quenching errors, at least on the 
B-meson decays constants, may be reduced by normalizing these constants with /^j^. The reduction of quenching 
effects is about 7% and the remaining effects become consistent with zero. To the extent that quenching errors on the 



B parameters are negligible, all of these considerations carry over to the quantities in Table IX which are proportional 
to decay constants. It should be remembered, however, that the real world has a third dynamical light quark and the 
effects discussed above may be amplified. 

A thorough estimate of quenching effects for all of the quantities that we calculate would require a dedicated 
unquenched simulation, which is beyond the scope of this paper. Therefore, we do not attempt to quantify these 
effects specifically. Nevertheless, as we mentioned in the two previous sections, the uncertainties in the lattice scale 
and in the strange quark mass are, at least in part, quenching effects. 

VIII. FINAL RESULTS AND DISCUSSION 

Because of the excellent consistency of the results for Bb^, Bb^, BbJBb^ and JbJ fsa &t the two lattice spacings 
and because we cannot perform a continuum extrapolation with only two points, we take the results from the finer 
lattice at /3 = 6.2, which should have smaller discretization and matching uncertainties, to be our best estimates. 

For the S'C/(3)-breaking ratio r^^, we also have a choice between the "direct" and "indirect" approaches described 
in the Introduction. While both methods give results that are compatible at the two values of the lattice spacing, the 
"direct" method leads to larger variations with lattice spacing and significantly larger statistical errors. Furthermore, 
as discussed in Sections^ and VI, the light- and heavy-quark-mass extrapolations are better behaved in the "indirect" 



method. Therefore, we take r^d obtained with the "indirect" method at = 6.2 as our best estimate for this quantity. 

A summary of our results for quantities directly relevant for B—B mixing is given in Eq. (^). These results are 
compatible with previous calculations of some or all of these quantities, which were performed using less improved 
relativistic fermion actions jl^, as well as with the recent calculation of Ref. |2J|, which makes use of a non- 
perturbatively, O (a)-improved, Sheikholeslami-Wohlert action. While the decay constants in Ref. |Q are non- 
perturbatively improved, B parameters and four-quark matrix elements are not. Thus, for those quantities, the 



discretization accuracy of that calculation is formally the same as ours. Moreover, the authors of Ref. |24| do not 
investigate cutoff dependence, as we do here with our two lattice spacings. Our results are also consistent, once 
systematic errors are taken into account, with the NRQCD results of Refs. |25,^, whose B parameters are 7 ^ 10% 
smaller than ours. 

For the decay constants and quantities proportional to them, the situation is less favorable than for B parameters and 
5t/(3)-breaking ratios. We do observe a two-statistical-standard-deviation dependence on lattice spacing, indicating 
that discretization errors are more important here. We quantify these discretization effects, as described in Section 



VII A. The corresponding uncertainty at /3 = 6.0 is large enough to bridge the gap between the results for i3- meson 
decay constants at the two lattice spacings. For D-meson decay constants, agreement requires that one also take into 
account the statistical error on the /3 = 6.2 results. So we take as our best estimates the results from the finer lattice, 
which in principle have smaller discretization and matching uncertainties, confident that our errors are a reasonable 
estimate of the uncertainty associated with this fixed lattice spacing calculation. In the future, though, when fo^ is 
accurately measured experimentally, it will be advantageous to consider the values of these quantities in units of this 
decay constant. With this normalization, the discrepencies between the results a.t (3 — 6. 2 an d 6.0, as well as the size 



of systematic (and statistical) errors, are significantly reduced, as can be seen in Tables VII and IX. 

Our results for the decay constants are summarized in Eq. (^. They are compatible with other recent calculations 
in the quenched approximation p2|- ^ , |3^ ,^, as reviewed in Refs. |4^, p6| , |47| ] . A very recent, non-perturbatively O (a)- 
improved calculation by the UKQCD Collaboration js^ yields Jb^ and Jb which are over two statistical standard 
deviations higher than our results. The authors of Ref. use the scale tq ||7^ to set the inverse lattice spacing, 
a~^, which gives an at the top of our range. They further use the preliminary non-perturbative values of 6a 
obtained in |7^,Q and the non-perturbative values of ca obtained in |79] . R Agreement with our results is recovered. 



Almost: the authors actually use the perturbative value of 6a. 
^®The authors of Ref. have since finalized their determinations of Ba, as well as those of ca The authors of Ref. 

are currently reviewing their analysis of decay constants to incorporate, amongst other considerations, this new information [Kl| 
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nonetheless, when systematic errors are considered. For _D-meson decay constants, their results agree with ours within 
statistical errors. 

All of our results are obtained in the quenched approximation. Some of the error resulting from this approximation 
is accounted for by the variations that the uncertainties in the lattice spacing and the strange quark mass induce, 
since these uncertainties are, at least in part, quenching effects. However, a thorough estimate of quenching effects 
for all the quantities that we consider here would require a dedicated unquenched simulation, which is beyond the 
scope of this paper. 



IX. CONCLUSION 



We have reported on high statistics, quenched lattice calculations of matrix elements relevant for neutral meson 
mixing, and leptonic decays of B and D mesons. We use mean-field-improved Sheikholeslami-Wohlert actions to 
describe quarks and work at two values of the lattice spacing. We have performed an extensive study of systematic 
errors and we believe that our final results, presented in Eqs. (||), (|^) and (p^), carry errors that reflect conservatively 
the uncertainty associated with our fixed lattice spacing calculations. 

Our results for neutral _B-meson mixing are compatible with the results of other calculations of some or all of the 
quantities we consider [T^,[2^-p^, as well as with the world averages of ||l7|,|4^-^. The same is true of our results for 
the decay constants, which are compatible with other modern, quenched determinations |^^-|36|, ^,|39|] and the world 
averages of ^^-Q. 

Finally, it should be emphasized that all of these results are obtained in the quenched approximation. They include 
a quenching error only to the extent that the variations in the lattice spacing and strange quark mass that we account 
for are quenching effects. It is worth noting that when dimensionless quantities are considered, such as B parameters, 
S'?7(3)-breaking ratios or quantities normalized by /d^, a number of systemaitc uncertainties, including possibly those 
associated with quenching, partially cancel. Nevertheless, a thorough quantification of quenching effects for neutral 
_B-meson mixing matrix elements would require a dedicated unquenched simulation, which is beyond the scope of 
this paper. It is important, however, that such a study be undertaken. The pioneering Uf ~ 2 studies of decay 
constants p8|-^ are a first step in this direction. 
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TABLE L Simulation parameters. Csw is the clover coefficient. 



13 6.2 6.0 

Csw 1.442 1.479 

size 243 X 48 16^ x 48 

# of configurations 188 498 



TABLE IL Hopping parameters, k^, of the light {tj} = q) and heavy (1/) = Q) quarks used in the simulations. In brackets, we 
give the masses, m^, of light-light pseudoscalar mesons composed of a degenerate quark and antiquark with hopping parameter 
Kq, obtained as described in Section We also provide the masses, Mp, of heavy-light pseudoscalar mesons composed of a 
heavy quark with hopping parameter kq and a massless antiquark, obtained as described in Section IV and Section Error 
bars are statistical only. The scales used to translate these masses into MeV are = 2.73 GeV at /3 = 6.2 and — 2.00 GeV 
at /3 = 6.0, as obtained in Section ^ 



~^ Kq [fUp in MeV] kq [Mp in MeV] 

6.2 0.13640 [831(4)], 0.13710 [608(4)], 0.120 [2238(14)], 0.123 [2006(13)], 

0.13745 [466(6)] 0.126[1757(11)], 0.129 [1488(10)] 

0.132[1186(9)] 

6.0 0.13700 [827(1)], 0.13810 [587(3)], 0.114 [2183(6)], 0.118 [1971(5)], 

0.13856 [454(3)] 0.122 [1746(4)], 0.126 [1503(4)], 
0.130 [1234(3)] 
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TABLE III. The coupling 0^3(9) at different values of q as obtained in the simulations at /? = 6.0 and 6.2. 



/3 


aMs(l/a) 


a!Ss(2/a) 


Q;Ms(7r/a) 


6.2 


0.1730 


0.1402 


0.1250 


6.0 


0.1921 


0.1522 


0.1343 



TABLE IV. Fit ranges for the correlation finictious and ratios used in this work. 



(3 6.2 6.0 

C^'^" /CfpCfp M<t^< 38 33 < < 37 
10 < < 14 11 < < 15 
C^'^"' /C^PpC^PP M<t^< 38 ■ii<t^< 37 

10 < < 14 11 < < 15 

C^PA{ix)/C^Pp{t^) 15 < < 23 13 < < 23 
C^'^ 13 < < 23 ll<t^< 23 

Cpp 13 < < 23 11 < < 23 

g^^p^ 10 < < 20 5 < < 23 

C^ajtS/C^p itx) 15 < < 23 11 < < 23 



TABLE V. Lattice spacings, critical hopping parameters and bare, strange quark masses, obtained as detailed in the text. 
Errors are statistical. 



f3 


6.2 6.0 


a-i (GeV) 


2.73 tf 2.00 tt 




0.13792(2) 0.13921(1) 


arhs 


0.0280(19) 0.0401(17) 
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TABLE VI. Results for the fit parameters for the heavy-quark-mass dependence of the various quantities studied. Cx = 
indicates a linear fit and Cx = Bx = 0, a fit to a constant. 



13 Ax Bx Cx 




6.2 0.94(2) -0.08(1) 
6.0 0.95(2) -0.10(2) 




6.2 0.96(4) -0.10(2) 
6.0 0.93(5) -0.12(3) 


Bp. 1 Bp^ 


6.2 0.97(3) 0.02(1) 
6.0 1.01(4) 0.03(2) 




6.2 0.010(1) -0.29(1) 
6.0 0.058(4) -0.85(2) 0.22(1) 




6.2 0.006(1) -0.29(2) 
6.0 0.041(6) -0.83(6) 0.20(3) 


{{p,\or'-''\p.)i{Pd\or 

x{MpjMp^) 


='^\Pd)) 6.2 1.57(26) 0.01(6) 
6.0 1.28(14) 0.11(8) 


^fiMp,) 


6.2 0.10(1) -0.32(4) 0.04(1) 
6.0 0.19(1) -0.53(2) 0.11(1) 




6.2 0.09(2) -0.30(8) 0.04(2) 
6.0 0.17(1) -0.51(4) 0.11(2) 


(/Pe//Pj 

Xy^MpjMp 


6.2 1.16(2) 
6.0 1.15(1) 



TABLE VIL Results at the two values of the lattice spacing. The first error on each quantity is statistical while the others 
correspond, respectively, to the variations in procedure described in the first five subsections in Section V^IL The errors enclosed 
in brackets reflect the variations considered in quantifying discretization errors. 



P 6.2 6.0 



fB 


[MeV] 




+8 
-5 


+11 

--0 


+6 
-3 


+0 

-6 


+16 +0 

-16 -0 


+9 

205 _g 


+14 +23 

-9 -0 


+9 

-4 


+0 +18 +0 
-11 -18 -0 




fB 

Id. 


0.71 tl 


" +2 
-1 


+3 ' 
-0 


+0 - 
-0 - 


fO - 
-2 


f3 +0 
-3 -2 


0.76 t 


3 
3 


" +2 +6 

-1 -0 


1 +° 

-0 


+0 +3 +0 

-4 -3 -2 


fB. 


[MeV] 204:^^2 


+10 

-6 


+13 
-0 


1 +^ 

-3 


+0 

-7 


+14 +6 

-14 -0 


233 +5 


-\ 


-17 +27 " 

-11 -0 


+10 
-4 


+0 +16 +6 
-13 -16 -0 




1b. 
Id. 


0.82 +3 


" +2 

-1 


+4 " 

-0 


+0 - 
-0 - 


fO - 
-3 


f2 +0 
-2 -0 


0.87 t 


1 
1 


" +3 +7 

-1 -0 


+0 

-0 


+0 +2 +0 

-5 -2 -0 




fB. 
fB 


1.15 


' +0 

-0 


+0 " 

-0 


+0 - 
-0 - 


fO - 
-0 


f3 +3 
-2 -0 


L14 + 


1 
1 


' +0 +0 

-0 -0 


+0 

-0 


+0 +2 +3 

-0 -2 -0 


fD 


[MeV] 


210 tl' 


+5 
~4 


+4 " 
-0 


+8 
-3 


+0 
-0 


+14 +0 
-15 -0 


+4 

228 




+7 +8" 

-7 -0 


+10 

-4 


+0 +14 +0 
-0 -15 -0 


fD. 


[MeV] 


236 +g 


+6 

-5 


+4 T 

-0 


+9 - 
-3 - 


^0 - 
-0 - 


fl2 +7 
-12 -0 


254 +^ 




+9 +9" 

-8 -0 


+11 

-5 


+0 +11 +7 
-0 -12 -0 




fn, 

fD 


1.13 


" +0 

-0 


+0 " 

-0 


+0 - 
-0 - 


fO - 
-0 


f3 +3 

-2 -0 


T11 + 


1 
1 


" +0 +0 

-0 -0 


+0 

-0 


+0 +2 +3 

-0 -2 -0 
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TABLE VIII. Results at the two values of the lattice spacing (continued), rfj""" = (B^\0^'^=^\B^)/{Ba\0^'^=^\Ba). An 
error of x indicates that the variation has not been explicitly performed, but that it is believed to be small. An error of y 
means that the variation has not been explicitely performed. 



6.2 



6.0 



Rn (Ma) 


91 '^^ 


+x 

— X 


+0 
-0 


+4 
-0 


+x 

— X 


+0 
-0 


+0 
-0 


0.89 


+4 
-4 


+x 

— X 


+0 
-0 


+3 
-0 


+0 
-0 


+0 
-0 


+0 
-0 


BbAMb) 


+2 

0.90 _2 


+x 

— X 


+0 " 

-0 


+3 

-0 


+x 

— X 


+0 

-0 


+0 

-0 


0.91 


+2 
-2 


+x 

— X 


+0 ' 

-0 


+3 

-0 


+0 

-0 


+1 

-1 


+0 

-0 


Bd{Md) 


0.82 +3 


+X 
— X 


+0 " 

-0 


+3 
-0 


+x 
—x 


+1 
-1 


+0 
-0 


0.81 


+3 
-3 


+x 

— X 


+0 ' 

-0 


+2 
-0 


+0 
-0 


+1 
-1 


+0 
-0 




0.98 tl 


+x 

— X 


+0 ' 

-0 


+0 
-2 


+x 

— X 


+0 
-0 


+0 
-0 


1.02 


+3 
-3 


+x 

— X 


+0 ' 

-0 


+0 
-0 


+0 
-0 


+0 
-0 


+0 
-0 


^direct 


1.61 


+X 
—X 


+o" 

-0 


+10 
-0 


+y 

— X 


+9 
-7 


+12 
-0 


1.36 


+11 
-12 


' +x 

—X 


+0 
-0 


+1 

-0 


+6 
-0 


+8 
-0 


+7 

-0 



TABLE IX. Results at the two values of the lattice spacing (continued). r^'J''"^'"'* = ( Ijf^ 



6.2 



6.0 



fB,^: [MeV] 210 



+20 
-21 



0.89 



+7 



fB^^I^JMeW] 241 + 

fSa . / Rnlo 

7^V 



+14 
14 



1.02 



+4 
-4 



fD\lB^^° [MeV] 222 



+10 
-10 



+10 +13 
-6 -0 

+2 +4 
-1 -0 

+12 +15 
-7 -0 

+3 +5 
-1 -0 

+5 +5 
-4 -0 



+12 +0 +19 +0 
-0 -7 -18 -0 

+2 +0 +3 +0 
-0 -3 -4 -2 

+13 +0 +17 +7 

-1 -9 -16 -0 

+2 +0 +2 +0 
-0 -4 -2 -0 

+13 +0 +14 +0 
-0 -0 -14 -0 



240 



+12 
-11 



0.95 



+5 



277 



+7 
-6 



1.09 



+2 
-2 



240 



+7 



+16 +27 
-11 -0 

+3 +7 
-1 -0 

+20 +32 
-13 -0 

+4 +9 
-2 -0 

+8 +8 
-7 -0 



+15 +0 +21 +0 
-5 -13 -21 -0 

+2 +0 +4 +0 
-0 -5 -4 -2 

+18 +0 +18 +8 

-5 -15 -18 -0 

+2 +0 +2 +0 
-0 -6 -2 -0 

+13 +0 +14 +0 
-3 -0 -14 -0 



Sds V 


0.94 


+3 
-3 


' +0 

-0 


+1 ' 

-0 


+2 

-0 


+0 

-0 


+1 

-1 


+0 

-3 


0.94 


+2 
-2 


' +0 

-0 


+0 ' 

-0 


+1 

-0 


+0 

-0 


+1 

-1 


+0 

-2 




1.15 


+2 
-2 


' +0 
-0 


+o' 

-0 


+0 

-1 


+0 
-0 


+2 
-2 


+3 
-0 


1.16 


+2 
-2 


' +0 

-0 


+0 ■ 

-0 


+0 
-0 


+0 
-0 


+3 
-2 


+3 
-0 


^indirect 


1.38 


+6 
-5 


" +1 
-0 


+0 ' 

-0 


+0 
-3 


+0 

-1 


+6 
-5 


+8 
-0 


1.39 


+5 
-5 


' +1 

-0 


+0 ' 

-0 


+0 
-0 


+0 
-0 


+6 
-5 


+8 
-0 
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matched at 2/a 



/Cq = 0.1200 
ic =0.1364 
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20 



25 
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FIG. 1. Plateau for the B-parameter of the operator O^^ renormahzed at fj, — 2/ a in the NDR-MS scheme, for /3 = 6.2 
(i.e. Eq. @ times 3/{8Z\)). The correlation function is obtained for W < ty < 14 and 10 < T - < 14. Points with the 
same T — tx + ty are shifted for clarity. 




0.02 0.04 0.06 0.08 0.1 0.05 0.1 

a(mj + mg) a(mj + mg) 

FIG. 2. Fits of the squared, light-light pseudoscalar meson masses versus light-quark mass to the PCAC relation of Eq. 
@, at /3 = 6.2 and 6.0. 
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K:m =0 
p 






'f^-1.19(2) _ 


(amj^)2 = 


=0.0609(25) ' 
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p/ \ ----p/ 

FIG. 3. Interpolation of afp, according to Eq. (|3l|), which is used to determine the scale a~^, a,t f3 — 6.2 and 6.0, from fx 
and niK, as described in the text. 
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FIG. 4. Light-quark-mass dependence of the heavy-light pseudoscalar meson mass, Mp^ , and extrapolation and interpolation 
to vanishing quark mass and strange-quark mass at /? = 6.2 and 6.0. 
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arriq am^ 

FIG. 5. Light-quark-mass dependence of the heavy-Ught B-parameter, Bp^[5 GeV), and extrapolation and interpolation to 
vanishing quark mass and strange-quark mass at /3 = 6.2 and 6.0. Points with the same arhq are shifted for clarity. 
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FIG. 11. Lattice results for X(Mp,) = $^(Mp,) versus l/(aMp,) at /? = 6.2 and 6.0 and for g = s and I. The solid curves 
are quadratic fits of the heavy-quark- mass dependence given in Eq. (5^). 




1 2 0.0 0.5 1.0 1.5 2.0 

l/(aMp_) l/(aMp_) 



FIG. 12. Lattice results for X(Mp, ) = Bp^/Bp, versus l/(aMp, ) at /3 = 6.2 and 6.0. The solid line is a fit to the linear 
part of the heavy-quark-mass dependence given in Eq. (^). 
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the constant part of the heavy-quark-mass dependence given in Eq. (pq). 
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